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Abstract. Recent calculations of particle acceleration in su- 
pernova remnants (SNRs) are used to estimate the associated 
7-ray production. For source spectra which are power-laws in 
momentum (or rigidity) the production efficiency of 7-rays 
with energy Ej > 100 MeV is shown to be about a factor 
2-3 lower than the value conventionally used for the interstel- 
lar medium and to depend only weakly on the spectral index 
of the power- law (in the range expected). Because the energy 
transferred to accelerated particles is rather tightly constrained 
by the total Galactic cosmic ray power, if SNRs are the main 
source of Galactic cosmic rays, this leads to an almost model- 
independent prediction of the SNR 7-ray luminosity in the 
band > 100 MeV. 

A detailed discussion of instrumental sensitivities and back- 
grounds shows that detection of SNRs in the Ej > 100 MeV 
band with, for example, the Energetic Gamma Ray Experiment 
Telescope (EGRET) will be difficult, but should not be impos- 
sible. However, and significantly, the prospects look much bet- 
ter in the TeV band accessible to modern imaging atmospheric 
Cherenkov telescopes. It should be possible to detect SNRs out 
to distances of 10 kpc if the region of the ISM into which they 
are expanding has a high enough density (n > 0.1 cm~^) so 
that their 7-ray luminosity is high enough. 

Finally, it is pointed out that existing and planned air- 
shower arrays can place important limits on the extension of 
the accelerated particle spectra in SNRs to energies above 
100 TeV. In conjunction with spectral measurements in the 
TeV region and detections or upper limits in the 100 MeV band 
this could provide a crucial test of current theories of particle 
acceleration in SNRs. 
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1. Introduction 



is mainly based on two facts; first, that a theoretical acceler- 
ation mechanism exists with certain attractive properties (dif- 
fusive shock acceleration: two recent reviews which emphasise 
the astrophysical applications are Blandford and Eichler, 1987; 
Berezhko and Krymsky, 1988), and secondly, that the galactic 
SNae (and thus the resultant SNRs) are almost the only known 
potential sources with the necessary amount of available en- 
ergy. Although plausible, these are indirect and unsatisfactory 
arguments. It is clearly desirable to seek direct evidence for 
accelerated nuclei in SNRs. The best hope would appear to 
be gamma-ray observations. Indeed some thirty years ago it 
was this, in the context of cosmic ray sources in general rather 
than SNRs in particular, which was at the heart of the ratio- 
nale for missions like SAS-II and COS-B. If the nuclear com- 
ponent of the cosmic rays is strongly enhanced inside SNRs, 
then through nuclear collisions leading to pion production and 
subsequent decay, 7-rays will inevitably be produced^. In this 
paper we discuss some estimates of the 7-ray visibility of SNRs 
in the energy range around 100 MeV typical of the EGRET in- 
strument on the Compton Observatory, the TeV energy range 
accessible to ground-based Cherenkov telescopes and the range 
above 10 TeV which can be covered by modern airshower ar- 
rays. 

Diffusive shock acceleration, as applied to cosmic ray pro- 
duction in SNRs has been studied by many authors (eg 
Blandford and Ostriker, 1980; Krymsky and Petukhov, 1981; 
Prischep and Ptuskin, 1981; Bogdan and Volk, 1983; Moraal 
and Axford, 1983; Lagage and Cesarsky, 1983; Jokipii and Ko, 
1987; Berezhko and Krymsky, 1988; Volk, Zank and Zank, 
1988; Dorfi, 1990; Kang and Jones, 1991). In general, it must 
be said that all calculations contain ad hoc parametrizations of 
various aspects of the physics, especially the injection of par- 
ticles into the acceleration process. Recent comparisons (Kang 
and Drury, 1992) show that the 'simplified models' (Drury, 
Markiewicz and Volk, 1989 [DMV]; Markiewicz, Drury and 
Volk 1990) are good approximations to the much more de- 
tailed hydrodynamical calculations, and in particular that the 



There is no clear observational proof that the nuclear compo- 
nent of cosmic rays is produced in supernova remnants (SNRs) 
although this is widely believed, at least for particle energies 
less than about 10^'' eV per nucleon. This general conviction 
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^ There will also be Bremsstrahlung and Inverse Compton 
components, produced by the electrons, both primary (directly 
accelerated) and secondary (from vr* decays), known from ra- 
dio synchrotron observations to be present in SNRs. By con- 
sidering only the gamma rays of hadronic origin we obtain a 
lower bound to the overall 7-ray emission of SNRs. 
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differences (in the worst cases factors of order two in the energy 
densities) arc loss than the uncertainties in the paramotriza- 
tions of the underlying physics. This suggests that they can be 
used to obtain estimates of the gamma-ray flux to be expected 
from nuclear interactions in a SNR. 

It is important to note that although the injection rate has 
to be specified in an essentially arbitrary fashion, the overall 
rate is rather tightly constrained by the obsorvationally in- 
ferred energy requirements. When allowance is made for the 
energy-dependent escape of particles above a few GeV (Swordy 
et al, 1990) the traditional estimates of production efficiency 
arc raised from a few percent to something above 10% per su- 
pernova (Volk et al, 1985; DMV) if the bulk of the cosmic ray 
nuclei are indeed accelerated in SNRs. 



2. Expected 7-ray fluxes 

In the simplified models the SNR is divided into three regions: 

an interior region filled with hot gas and accelerated particles 
but very little mass, an immediate post shock region where 
most of the matter is concentrated, and a shock precursor 
region where the accelerated particles diffusing ahead of the 
shock compress the ambient medium. In each of these regions 
we know the energy density of the accelerated cosmic ray par- 
ticles and the mass density. The production rate of 7-rays per 
unit volume can be written as 



(1) 



where n is the number density of the gas, Ec is the CR energy 
density and is the ">-ray emissivity normalized to the CR 
energy density, q-y = E^/Ec- For the galactic CR in the diffuse 
interstellar medium (ISM), outside sources, the 7-ray emis- 
sivity above 100 McV due to nuclear collisions is about E-yi^: 



(e.g. Stecker, 
The corre- 



100 MeV) ^ (1 - 1.5) X 10"^-'s"'(H-atorn)" 
1971; Dernier, 1986) and Ec ^ 10"^^ erg c: 
sponding q^{> 100 MeV) « 10"^^ s"^erg"^cm^(H-atom)"^ in 
the ISM can also be used for rough estimates of the 7-ray lu- 
minosity of SNRs (see below). 

In fact, because the spectrum in the source is expected to 
be harder than the observed cosmic ray spectrum (which we 
know, from the energy dependence of the relative abundance 
of spallation products, to have been softened by propagation 
effects; JuUusson et al, 1972) this is almost certainly an overes- 
timate of the production rate in SNRs. The spatially averaged 
spectrum in the source is probably of the form /(p) (x p^" 
with 4 < a < 4.3, whereas the observed cosmic ray spectrum 
has a ~ 4.7. Gamma rays above 100 MeV are produced by 
particles of energy above about 1 GeV per nucleon and the 
cross-sections are not very energy-dependent; thus the gamma 
production rate may be taken to be proportional to the num- 
ber density of relativistic particles. For a given energy density, 
the harder the spectrum the lower will be the total number 
density of paxticles. On the other hand, if the spectrum ex- 
tends to energies below about 1 GeV per nucleon, a harder 
spectrum will have fewer particles which contribute to the to- 
tal energy density but do not generate observable 7-rays. To 
quantify these competing effects let us, following DMV, con- 
sider a simple power-law spectrum of the form 



where / is the (isotropic) phase space density and p is momen- 
tum with a lower cut-off at pinj and an upper cut-off at Pmax. 
Then the number density of relativistic particles with p > 2mc 
is easily calculated to be 



Nr 



1 



a — 3 



47r/o(2mc)^ 



(Pmax \ ^ 
mc ) 



(3) 



and for a > 4 and pmax S> mc the dependence on the upper 
cut-off can be neglected. On dividing by the corresponding 
energy density (as calculated in DMV) we obtain 



Ec 2mc2 



2mc V 2mc / _ 
- (pinj/2mc)5-« 



2(a - 3)mc2 

(pinax/'2/;i(:)^ 



+ - 
1 



5 — a 



- 1 



(4) 



4 < a < 5, 



4mc^ 



, / 2mc 
In 



1 - 



2mc 



It is easy to see that for typical values of pinj/mc ~ 10 ^ — 
10"^ and pmax/mc « 10" - 10® the ratio Nr/Ec is rather 
independent of a even in the most extreme case of a hard source 
spectrum proportional to p~* and an ambient spectrum of the 
form p-" '^. This semi-quantitative argument is confirmed by 
rmmcrical calculations which show that for the spectrum of 
accelerated particles represented by Eq. (2) the production rate 
is almost independent of the spectral index a. In fact, for 4.1 < 
a < 4.7 the production rate may be considered constant within 
20% (see Table 1) of 



q^{> 100 MeV) « 0.5 x 10"^^ s"^erg"^cm^(H-atom) 



(5) 



which is a factor of three less than the corresponding value 
for a Galactic CR flux of the form Jj,{E) (x {E + mc^)-'^ ''^ 
where E is kinetic energy. This particular form for the spec- 
trum gives a very high (close indeed to the maximum possible) 
efficiency for the production of > 100 MeV 7-rays (at fixed CR 
energy density) since it contains negligible energy density in 
nonrelativistic particles (useless for gamma production) and, 
as mentioned above, at the same time heis a low energy density 
for particles with energy much above a few GeV. As the 7r°- 
production cross-section becomes constant above a few GeV 
the maximum value is reached in the case of a monoenergetic 
beam of protons with energy 1-2 GeV. The results of numeri- 
cal calculations of qy{> 100 MeV) for different spectral indices 
a are presented in Table 1. 

At energies E ^ 100 MeV the production spectrum of 7- 
rays reproduces the spectral shape of the distribution of parent 
protons. The emissivity, for power- law spectra of relativistic 
protons can be presented in a simple analytic form (Appendix 
A). In the case of a source spectrum of the form given in Eq. (2), 
the production rate, normalized to the energy density of the 
accelerated particles is 



/(p)=/o(p/2mc)- 



(2) 97 (> E) 



(6) 
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Table 1. The production rates of 7-rays for different spectral 

indices a of the parent cosmic ray distributions, f(p) oc p^" ■ 
The contribution of nuclei other than H in both the CR and 
the target matter is assumed to be the same as in the ISM, 
ie the pure proton contribution is multiplied by a factor 1.5 
(cf Dermer, 1986). In calculating the low energy part of the 
production rate we used the formalism developed by Stecker 
(1971). The units are s"^erg"icm^(H-atom)"^ 



a g.y(>100MeV) g^(>lTeV) 



4.1 


0.46 X 10" 


13 


1.02 X 10 


-17 


4.2 


0.58 X 10" 


-13 


4.9 X 10" 


18 


4.3 


0.61 X 10" 


-13 


2.1 X 10" 


18 


4.4 


0.57 X 10" 


-13 


8.1 X 10" 


19 


4.5 


0.51 X 10" 


-13 


3.0 X 10" 


19 


4.6 


0.44 X 10" 


-13 


1.0 X 10" 


19 


4.7 


0.39 X 10" 


■13 


3.7 X 10" 


20 



The total gamma-ray luminosity is given by j q-yuEc d'^r 
which, in the context of the simplified models, can be written 
as q-y {MiEci + M2EC2) where Mi is the total mass in the 
precursor region, M2 that in the immediate post-shock region 
and Eci,2 are the corresponding cosmic ray energy densities. 
Because of diffusion across the shock the two energy densities 
arc in fact equal, Eci = Ec2- The common value is also not 
very different from Ec:i, the energy density in the interior of 
the remnant, for two reasons. First, there is diffusive coupling 
between the acceleration region around the shock and the in- 
terior of the remnant. Secondly, if the acceleration is efficient, 
the cosmic rays provide a substantial, if not the dominant, part 
of the interior pressure and the interior of the remnant has, 
for dynamical reasons, to be in rough pressure equilibrium. 
It follows that, to order of magnitude, the cosmic ray energy 
density throughout the remnant and in the shock precursor is 
Eci — Ec2 ~ Ec:i ~ SOEsn/^ttR^ where 6 is the fraction of 
the total supernova explosion energy, Esn, converted to cos- 
mic ray energy and R is the remnant radius. Thus for the 7-ray 
luminosity we can estimate 



where q'y{> 1 TeV) is also tabulated in Table 1 for various 
spectral indices a. 

The simplified models do not as yet include any cooling 
of the gas (in contrast to the hydrodynamic models of Dorfi, 
1991a,b), however we do not believe that this is important for 
the estimates presented here. The cosmic ray energy density, 
because of efficient diffusion, is rather uniformly distributed 
inside the remnants so that clumping of the gas as a result 
of cooling is not likely significantly to change the results. Of 
course cooling does change the shock strength, and therefore 
the acceleration efficiency, in the late stages of SNR evolution. 

To discuss the visibility of SNR as 7-ray sources two quan- 
tities seem most interesting. The first is the total luminosity 
of the remnant in photons per second, the second the surface 
brightness along a line through the centre of the remnant in 
photons cm"^ s"'^ sr"^. While the first measure is appropriate 
to discussions of the detectability of unresolved sources, the 
second is more appropriate for spatially resolved sources. Fig- 
ure 1 shows both measures for a rather standard SNR model 
together with the radii of the shock, the thickness of the precur- 
sor and the partition of the explosion energy between thermal 
and nonthermal particles. The dotted vertical line indicates the 
point where the post-shock gas temperature drops below 10® K 
and cooling becomes potentially important (although as noted 
we have not included cooling in the present calculations). For 
this calculation the ejccta mass was set to 1 Mq and the am- 
bient density to 0.1 cm"'\ Alfven heating in the precursor and 
diffusive coupling between the interior of the remnant and the 
precursor were included with parameters an = 1 and ac = 10; 
the injection was parameterized by specifying that a constant 
fraction e — 10"^ of the kinetic energy dissipated in the shock 
was given to energetic particles capable of further acceleration. 
The physical processes represented by the parameters e, aH 
and ac are discussed in detail in DMV. Here we note simply 
that the effect quantified by ac is clearly required on theoret- 
ical grounds and has been shown to be a significant effect in 
recent numerical studies (eg Dorfi, 1991a,b). The importance 
of Alfvcnic effects is the subject of a recent study by .Jones 
(1993). It is worth noting that both effects act to reduce the 
amount of acceleration and thus tend to decrease the 7-ray 
production. 



= g^(Mi + M2) ~ eq.,EsNn 

« 10^®6I(Ssn/10^^ erg)(n/l cm"^) phs"^ 

where ri is the ambient density. Obviously the exact value of 
6 depends on the details of the model, in particular the in- 
jection parametrization used, but clearly < 1 and, as dis- 
cussed above, the observations appear to require > 0.1 at 
the point when acceleration stops and the produced cosmic 
rays diffuse out into the ISM. In fact for reasonably plausible 
injection models 6 is roughly constant throughout the Scdov 
phase with only a moderate dependence on external parame- 
ters such as the ambient density (Markiewicz et al, 1990). At 
earlier times, during the free-expansion phase, the amount of 
energy processed through shocks is small (almost all of the ex- 
plosion energy is in the form of kinetic energy of the ejecta) 
and 61 < 1. 

The calculated results agree well with this simple picture; 
the total luminosity is low during the free-expansion phase, but 
is roughly constant during the Scdov phase. The Scdov phase 
(or better, Sedov-like, because the models are not exact sim- 
ilarity solutions) starts when the amount of swept up matter 
is roughly equal to the ejecta mass. Typically this occurs at 
a radius 0(10 pc) and the SNR spends most of its useful life 
(until the outer shock becomes weak) in this phase. Cooling 
of the post-shock gas is important when the post-shock tem- 
perature drops below about 10® K, but is not as catastrophic 
for these models as for pure gas models because the CR do 
not cool and continue to provide an internal pressure in the 
remnant. As long as the outer bounding shock is strong and 
active in accelerating particles the diffusion coefficient for cos- 
mic rays near the shock front is very small (^ RR except at the 
highest energies). Thus the accelerated particles are confined 
to the remnant and only diffuse out into the ISM when the 
shock, either because of expansion or cooling, is no longer able 
to maintain the strong levels of magnetic disturbance which 
lead to diffusion coefficients very much smaller than those in 
the ambient ISM. 

If the SNR is a distance d away, this translates into a flux 
at the earth of 
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Fig. 1. The evolution, as a function of age in years, of the total 7-ray luminosity (in photons with energy above 100 MeV per second) and 
the central surface brightness (in photons cm~^sr~^s~^) together with the shock radius and the width of the precursor region (in pc) and 
the partition of the explosion energy between accelerated particles [Eq) and thermal gas {Eq) for a SNR expanding into a medium of 
density 0.1 cm~^. 



F(> 100 MeV) w 4 4 X 10"^6l [ -^ElS] ( ^ proportional to E or a = 4.1, the flux in the high energy 

^~ ' ''^ ' \ IQSi erg y \1 kpc J ^^g^ 7-ray region is 

Vlcm-sy 



-2 -1 
cm s 



1051 erg 

7-rays with energies above 100 MeV. Similar estimates have 2 ^ (9) 

been given by Dorfi (1991a,b). For the same SNR parameters, 



and assuming a differential energy spectrum inside the remnant V 1 kpc 



(lcm-3) 
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Fig. 2. As Fig. 1, but for an external density of n = 5 cm ^. 

The surface brightness is highest approximately at the tran- 
sition from the free-expansion phase to the Sedov phase. For 

fixed ejecta mass the radius at which this occurs scales as n~^^^ 
so that, combining this with the earlier estimate of the total 
luminosity, we deduce that the maximum surface brightness 
should scale with the external density as n"'^'^. 

In Fig. 2 results are shown for an explosion in a higher 
density medium with n = 5 which confirm the above estimates 
(note that the neglect of cooling is more serious in this case). 
Although SNRs in denser regions of the ISM are stronger 7-ray 
sources than those in low density regions, they are probably not 



as easy to detect because of the problems associated with back- 
ground subtraction and source confusion. Indeed the COS-B 

measurements show that the diffuse E > 100 MeV 7-ray emis- 
sion at |6| < 5° is about 5x lO"** and 1.5x 10"'' phcm^^s^^sr"' 
from the inner and outer parts of the galactic plane (Bloemen, 
1987). Thus the expected flux of 7-rays within 1° (the typical 
size of the near-by SNRs) is about (1.5 — 5) x 10~^ making 
the detection of any source in the galactic plane with an in- 
tensity less than 10~^ phcm~^s~'^ hard even for the EGRET 
experiment (see Fig. 4). 
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Fig. 3. As Fig. 1, but for an injection model which places most of the injection and acceleration at late times (see DMV for details). 



It is interesting to ask whether, by allowing the injection to 
be weaker in strong shocks (as may be the case) , it is possible 
to reduce the surface brightness to the point whore SNRs would 
be hard to detect. In DMV a model of Mach number depen- 
dent injection was discussed which has the property that most 
of the acceleration occurs towards the end of the Sedov phase. 
In Fig. 3 the results for this model {ai = 1 and e = 0) are 
presented for n = 0.1. As expected the final total luminosity 
(which is determined essentially by the explosion energy and 
the ambient density if the acceleration is efficient) is not signif- 



icantly affected, but the surface luminosity is reduced because 
of the larger remnant size. 

Before turning to the observational implications of these 
models we note that they do assume expansion of the SNR into 
a uniform medium. While this is probably a good assumption 
for SNae of type la, in other cases we would expect the ambient 
medium to be perturbed by the stellar wind and photon flux 
of the SN precursor. In detail this has not been explored, but 
qualitatively we expect that there will be rather little CR pro- 
duction at early times as the young SNR expands into a low- 
density interstellar bubble of radius < 10 pc (although what 
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acceleration there is may be to energies above 10^* eV; c/ Volk 
and Biormann, 1988). Roughly two cases may have to be dis- 
tinguished depending on the precursor star. If the SNR has 
swept up a wind mass which is larger or of the same amount 
as the SN ejecta before reaching the edge of the bubble, then 
the internal energy density and the pressure are roughly uni- 
form over the SNR and therefore 'diluted' in the large bubble 
volume. With this postsliock pressure the outer shock will now 
shock heat the ISM matter swept up and compressed by the 
wind (and expanding Hll-region), increasing the CR produc- 
tion and L~f. However 9 may not reach the value which would 
be possible if the ambient ISM had not been perturbed. In the 
second case, for smaller mass of the precursor star, the wind 
mass is considerable smaller than the ejecta mass. When the 
latter masses reach the edge of the bubble and shock heat the 
swept-up wind shell, the subsequent evolution should be very 
similar to that of a SNR expanding into a uniform medium. 
The main effect is thus probably to delay any significant CR 
production (and associated 7-ray production) until the SNR 
has reached the edge of the wind cavity. 



3. Discussion 



of the CR are much more frequent, it is possible to discrimi- 
nate, at least partially, between hadronic and electromagnetic 
showers using various characteristics of the shower secondaries. 
At TeV energies morphological analysis of the Cherenkov light 
image has been shown to be a powerful technique for reject- 
ing hadronic showers, and at -E > 100 TeV the muon content 
can be used. The sensitivities of modern ground-based detec- 
tor systems to cosmic 7-ray fluxes are mainly determined by 
the efficiency of hadron rejection, which has to be estimated 
by detailed Monte Carlo simulations. 

The observational possibilities for detecting 7-rays from 
SNRs are summarized in Fig. 4. At low 7-ray energies, E > 
100 MeV, the typical range of the diffuse COS-B fluxes from 
the galactic disc (|6| < 5°) are shown, bracketed from above 
by the average emission from the inner Galaxy and from below 
by the slightly harder average emission from the outer Galaxy 
(Bloemen, 1987). We plot the flux that would bo observed by 
an instrument with an acceptance of 1°, roughly corresponding 
to the expected angular size of SNRs. The three upper curves 
are calculated spectra from SNR sources using our accelera- 
tion models for three source spectra, /(p) oc with a = 4.1, 
4.2 and 4.3, and using the omissivities from Appendix A. The 
spectra are normalized such that 



These results suggest that it should be possible to detect 

> 100 MeV 7-rays from at least some of the relatively near- 
by (d < 1 kpc) SNRs, especially at high galactic latitudes whore 
the background confusion is less. It is interesting to note that 
claims have been made for a detection of the north polar spur 
(thought to be an old near-by remnant) in the COS-B and 
SAS-II data at flux levels similar to those discussed here, but 
the significance level is very low and the results depend cru- 
cially on the background subtraction model used (Lebrun and 
Paul, 1985; Bhat et al, 1985). 

The other, perhaps more promising, possibility to test the 
SNR origin of cosmic rays, is to look for very high energy 7 rays 
from these objects. If, as the models suggest, the accelerated 
particle spectra in SNRs are hard [a ~ 4) and extend to at least 
10'^'* eV, then the fluxes of TeV 7-rays expected are well above 
the sensitivities of modern Cherenkov telescopes, even for SNR 
located in low density regions of the ISM. Equation (9) shows 
that, for Si?sN ~ 10^'^ erg, the flux of TeV 7-rays oxpoctod 
from a near-by SNR (at d < 1 kpc) exceeds 10"^^ phcm^'^s^^ 
at n > 0.1 cm^'\ 

Because the diffuse CR in the ISM have a softer spectrum 
than that expected for SNR source spectra the diffuse Tr^-decay 
background should be much less of a problem at TeV energies 
and beyond than at 100 MeV. There is also a diffuse Galactic 
7-ray background resulting from IC scattering of the 2.7 K 
cosmic microwave background by the CR electrons. But at TeV 
energies this also appears to be relatively less important than 
at S « 100 MeV (soo Appendix B). 

In fact for ground- based detectors, which arc the only ones 
available at energies E > 0.1 TeV, the principal background 
flux is that of the cosmic ray electrons themselves because they 
produce atmospheric showers which are indistinguishable from 
7-ray induced showers of the same primary energy^. Although 
the background showers induced by the nucleonic component 



A = I 



^ Of course this is not a problem for spacecraft instruments 
which can reject charged particle events by using an anti- 
coincidence shield. 



1Q51 g]-g / V 1 cm 



d 



1 kpc 



= 1 



(10) 



where the quantities are as defined in section 2. For 1 the 
plotted source spectra should be scaled by a factor A. Remark- 
ably these curves roughly coincide at 7-ray energies < 1 GeV, 
ic in the EGRET range. 

The diffuse 7-ray emission, again within 1°, for a total hy- 
drogen column density (both atomic and molecular) (A'h) = 
(iVH2) + (-'VHi) = 10^^ cm~^, a value typical of the outer Galaxy, 
and a cosmic ray confinement volume of size L = 10 kpc, are 
given by the two lower curves. The curve labeled 7r° corre- 
sponds to a diffuse Galactic cosmic ray spectrum of the form 
Jp{E) <x. i^E + rUpC^) ' (see section 2)^. 

The curve marked IC corresponds to the sum of the diffuse 
IC gamma ray emissions as discussed in Appendix B and above. 
Finally the background equivalence of the high energy diffuse 
interstellar electrons is given by the dash-dotted line, based on 
the observations of Nishimura et al (1980, 1990; see Appendix 
B). 

The four hatched curve segments above 0.1 TeV correspond 
to the 5cr sensitivities of different ground-based detection sys- 
tems for one year of operation (w 100 h and w 2000 h source 
exposure for a single imaging telescope and charged particle 
detector arrays respectively). The sensitivity (defined as the 
minimum flux of a 7-ray source detectable at a level of confi- 
dence of m standard deviations) is 



-,(inin) 



(> E): 



Pcit.S'cQ f Jcn{> E)AnY^' 
P 9 



(11) 



While the spectrum above 10 GcV is well determined, at 
lower energies, where solar modulation is dominant, we do not 
really know the form of the demodulated ISM spectrum. This 
particular form is compatible with the observations, and has 
been extensively used by many authors in the past (eg Stecker, 
1971), but otherwise has no special significance. 
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Fig. 4. The predicted 7-ray fluxes from a SNR together with the various Galactic backgrounds (for a 1° source) and current instrumental 
sensitivities; see text for details. 
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where Jcr(> E) is the integral cosmic ray flux, AQ ~ nAtf)'^ 
where A(f> is the angular resolution of the detector, S-y and Scr 
are the effective detection areas for photons and cosmic rays 
respectively (in general the two are not equal), PcR^-y is the 
rejection factor for background events (the probability that a 
CR event will be misclassified as a 7-ray event) and P-,^^ is the 
probability of correct classification of 7-ray induced showers. 
The angular and lateral distribution of the Clierenkov radiation 
as well as the muon content are characteristics which allow 
a rather efficient separation of hadronic and electromagnetic 
showers, PcR^-y -C 1 and Pj^j « 1 (eg Weekes, 1992). The 
high efficiency of the Cherenkov imaging technique has been 
demonstrated recently by the Whipple collaboration (Vacanti 
et al, 1991). 

It should be noted that the detector sensitivities shown in 
Fig.4 correspond to point sources, namely to sources with an- 
gular size smaller than A(j), where A(j) ~ 0.1° corresponds to 
the resolution of modern imaging Cherenkov telescopes. For 
extended sources the high resolution may not be of great help 
as far as nucleon rejection is concerned. In particular, in case 
of the imaging Cherenkov technique the 'orientational' param- 
eters of Cherenkov images become rather ineffective, and the 
separation of proton- and 7-induced showers should be done 
mainly through exploiting the 'shape' parameters . For a single 
telescope, the 'shape' parameters alone will provide a sensitiv- 
ity which is a factor 2 or 3 lower than that shown in Fig.4 (e.g. 
Reynolds et al, 1993). In fact, a system of imaging Cherenkov 
telescopes will allow a much more effective rejection of the 
CR background than a single telescope. The probability of de- 
tection of the 7-ray induced shower by a single telescope up 
to distances R ~ 100 m is close to 1, whereas the detection 
probability of CR protons near the effective energy threshold 
of 7-rays is rather low, namely < 0.1. Therefore the require- 
ment of detection of a shower by a system of several (> 3) 
telescopes,separated by 50 — 100 m from each other,practically 
does not reduce the number of detected 7-rays, but it "sup- 
presses" the CR background by factor of at least 10~^, i.e. at 
the effective energy threshold of 7-rays Sen < W~^Sj. More- 
over, the software analysis of the shape of the Cherenkov im- 
ages in different telescopes (different projections) will be of 
additional help for the improvement of the efficiency of 7/p 
separation (Aharonian et al.,1993). These two factors together 
allow us to expect that future systems of imaging Clierenkov 
telescopes will approach an almost 'background-free' detection 
(i.e. the number of detected CR is smaller than the number 
of detected 7-rays) down to 7-ray fluxes F-y ~ 10~^'^ph/cm^s, 
before the statistics of 7-rays becomes the limiting factor, even 
for extended sources. Nevertheless, the angular size of the 7-ray 
sources for a reasonable value of the aperture of the imaging 
camera (< 5°) should not exceed ~ 2°, since beyond 2° the 
efficiency of detection of 7-rays drops sharply. These results 
will be published elsewhere. 

As we will show below the radius of a SNR only exceeds 
10 pc after a minimum age of a few thousand years, provided 
we limit ourselves to ambient gas densities n > 0.1 cm~^. Ob- 
viously a SNR in a lower density environment would be very 
difficult to detect. Moreover for more favourable parameters, eg 
those characteristic of Tycho's SNR, the time when the source 
size reaches 10 pc is approximately 2 x lO** y. Note that a radius 
of 10 pc implies an angular size < 1° if the source is located 
at more than 1 kpc. On the other hand the 7-ray luminosity 
reaches a plateau level at ti < W* y, and continues to be at 



this level during several ii, depending on the parameters (see 
Fig. 1). Since F-y oc L^jS the best candidates for detection are 
SNRs with age of order 10'* years at a distance d r; 1 kpc. Of 
course for SNRs in dense environments the 7-ray fluxes may be 
considerable even for d > 1 kpc. In galaxies similar to ours the 
core collapse SNe form massive stars outnumber SN type la by 
almost an order of magnitude (Evans, van den Bergli, and Mc- 
Clure, 1989; Tammann, 1992), and such stars form typically as 
clusters in dense clouds. Given the influence of these stars on 
the ambient medium, it is nevertheless not clear how frequent 
SNRs in a n > 10cm~^ medium are. Still, the enhanced SNR 
rate in the molecular ring of the Galaxy at d ~ 4 kpc (eg Bloe- 
men et al., 1993) make this region in principle very attractive 
for the detection of distant SNRs. In the 100 McV region such 
detection appears practically quite difficult as becomes clear 
from the fact that the COS-B diffuse galactic fluxes within 
1° (a few 10~^ cm~'^s~* , see Fig. 4) approximately correspond 
to the sensitivity of the EGRET detector for discrete 7-ray 
sources in the galactic disc. 

As can be seen from Fig. 4, the situation is better at 
TeV energies, where 7-rays from SNRs may be detected up 
to d ~ 10 kpc provided n > Icrri^"'. (It is interesting to note 
that for SNRs with marginally detectable fluxes and angular 
sizes of order 1° a lowering of the energy threshold below 1 TeV 
is not necessarily the most efficient detection strategy because 
of the rising electron background.) Therefore future imaging 
Cherenkov telescopes (ICT) will be very useful in searching 
for high energy 7-rays from SNRs within wide ranges of ages 
and distances. The detection of near-by (d < 1 kpc) SNRs 
is complicated by the large angular size of the sources, but, 
as was mentioned above, possible with arrays of ICTs. The 
Whipple GRANITE detector (Akerloff et al, 1990) consist- 
ing of two « 10 m imaging telescopes has been built recently, 
and the HEGRA imaging telescope array (Aharonian et al, 
1991) consisting of five ICTs is now under construction at La 
Palma. The high angular resolution of ICT arrays will allow 
the 7-ray spectrum to be measured with rather high accuracy 
[AE / E ~ 20%) and thus give direct information about the 
spectrum of accelerated particles in the SNR. The weak depen- 
dence of the fluxes at 0.1 to 1 GeV on the shape of spectrum of 
accelerated particles (assuming a power-law form; see Fig. 4) 
means that measurements in this energy range arc in some 
ways less informative than those at TeV energies. However, for 
the same reason, EGRET measurements can provide a rather 
model-independent estimate of the normalization parameter A 
[Eq. (9)] and thus the total nonthermal power of SNRs. 

The important practical question of identifying candidate 
SNRs for observation can best be illustrated by an example. 
Let us consider Tycho's SN of 1572 which is at an estimated 
distance of d = 2.25 ± 0.25 kpc (e.g. Heavens, 1984). From 
the historical light curve the SN is commonly assumed to have 
been of type la, an accreting white dwarf system, for which our 
model assumptions of a uniform ambient medium are probably 
well justified (c/the discussion at the end of section 2). Accord- 
ing to the most recent radio observations (Tan and Gull, 1985) 
the expansion parameter ttt/R = 0.462 ± 0.024 where R{t) is 
the outer SNR radius and R{t) its time derivative at time t 
after the initial explosion. Thus the expansion parameter has 
decreased quite far towards its asymptotic value of 0.4 in the 
Sedov phase. This dynamical argument is probably more re- 
liable than present X-ray models which need to take account 
of non-equilibrium electron temperatures and non-equilibrium 
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ionisation states behind the outer shock as well as nonsolar 
chemical composition in the cjocta which arc part heated by 
the reverse shock. Thus Tycho's SNR is probably quite close to 
the Sedov phase. Different X-ray models give results ranging 
from n = 0.28 cm"^ and Esn = 7 x 10^° erg (Hamilton et al, 
1986) to n = 1.13 cm"^ and Ssn = 1.9 x 10®° erg (Smith et al, 
1988). A value of Esk/u = 0.17 x lO'''^ crgcrn^ with a quoted 
error of about 20% from the latter reference agrees reasonably 
well with the dynamical value of Ssn/w — 0.2 x 10*^ ergcm^ 
inferred by Heavens (1984) from radio and X-ray data and 
gives the value nSsN = 0.2 x lO^^ergcm"^. Thus assuming 
9 ~ 0.2 we have A ~ 0.01. If the Sedov phase is far from being 
reached, then we may have 6 -C 0.2 and A <C 10^^; however 
it would be hard to reconcile this with the amounts of shock 
heated gas seen in the X-ray band as well as the kinematical 
evidence from the expansion parameter. Alternatively, it may 
be that in this particular remnant particle acceleration is in- 
efficient and 6 <C 0.2; but then, if the Galactic SNRs are the 
sources of the cosmic rays, the remaining SNRs must be even 
stronger accelerators and we soon run into energy problems. 
Thus, on balance, the prediction of A m 0.01 for Tycho's SNR 
is a hypothesis worth testing. 

The spectrum inside such a young SNR is presumably very 
hard with a < 4.1. This is due to the fact that the shock is still 

very strong and wave heating is inefficient (DMV) so that even 
in a simple test particle approximation the spectrum is very 
hard (Volk, Zank and Zank, 1988). Nonlinear effects tend to 
flatten the high-energy end of the spectrum so that assuming a 
spectral slope of 4.1 for Tycho is not unreasonable. On looking 
at Fig. 4 we see that an EGRET detection of Tycho appears 
impossible, but that it should be visible at 1 TeV for imaging 
Cherenkov arrays and at 20 TeV for a 1 km^ AEROBICC array. 
The present radius of Tycho's SNR corresponds to about 4 arc 
minutes. 

This situation should persist for lO* years at which epoch 
Tycho's SNR will have an angular radius of about 0.25°. The 
fraction of the explosion energy converted to cosmic ray energy 
should remain at 6^ ~ 0.2 or higher while the spectral slope will 
begin to increase from a ~ 4.1 to a ~ 4.3 only towards the end 
of the evolution. Therefore a remnant of this type may be bet- 
ter observed at somewhat later times. It is worth noting that 
if Tycho were located at d = 1 kpc we would have ^ « 0.1 and 
detection would be much simpler. Obviously SNRs in denser 
environments are better candidates at the same distance and 
comparable evolutionary epochs. But we are dealing here with 
the statistics of a small number of nearby sources and each can- 
didate requires an individual evaluation before sensitive survey 
instruments become available. 

Most models of CR acceleration in SNRs predict upper cut- 
offs in the spectrum between 10^ and 10^ TeV. There will how- 
ever be some SNRs whose precursor stars were massive OB 
stars or Wolf-Rayct stars. These SNae explode into a medium 
where the magnetic field is that in the wind of the progenitor 
star and under these circumstances energies as high as 10^ or 
10* TeV may be possible (Volk and Biermann, 1991). If such 
special sources have power-law spectra a clear observational 
test of their existence is possible with the existing UMC (Utah, 
Michigan, Chicago) detector (Pick et al, 1991). The 5a sensi- 
tivity of this system above 100 TeV at present reaches about 
3 x 10~^* phcm~^s~^ in one year of operations (see Fig. 4). 
Therefore the efficiency of CR acceleration in SNR to energies 



significantly above 10^ TeV may soon be tested for all near-by 
SNRs using the all sky survey of the UMC detector. 

The new generation air shower detectors like CRT (the 
Cosmic Ray Tracking experiment; Heintze et al, 1989) MILA- 
GRO (Multiple Image Los Alamos Gamma Ray Observatory; 
Williams et al 1991) and the 1 km^ AIROBICC (AIR shower 
Observation By angle Integrating Cherenkov Counters; Lorenz, 
1992) could observe SNR sources above 10 TeV provided the 
source CR spectra extend to energies > 10^ TeV (note that 
the spectra in Fig. 4 have been drawn with no high energy 
cutoff). Thus even negative results, in the form of firm upper 
limits in the E > 10 TeV range, could be of great importance 
in determining the upper cutoff of the acceleration process in 
SNRs. 

Although this paper concentrates, for obvious reasons, on 
the 7-ray visibility, there is another possibility to detect high 
energy stable neutral particles from SNRs. From Table Al we 
expect (for hard source spectra, a « 2) a neutrino flux above 
1 TeV of 

F^^{> ITeV) w 10"^°^ neutrinos cm"^s"^ (12) 

Such fluxes may be marginally detectable for DUMAND (Deep 
Underwater Muon And Neutrino Detector; eg Okada, 1991) in 
the case of near-by SNRs. This is especially interesting because 
of the all sky survey capability of DUMAND. 

Besides SNRs in our own Galaxy, it may be possible to use 
the same combination of techniques to detect the integrated 
7-ray flux of some favourable near-by normal galaxies. Such 
galaxies do not have an active nucleus, but presumably con- 
tain the same types of cosmic ray sources as our Galaxy. Exam- 
ples might be the Magellanic Clouds and the near-by starburst 
galaxy M82. This will be discussed in a separate paper. 



4. Summary 

Recent nonlinear particle acceleration models for SNRs have 
been used to predict their 7-ray luminosities. A detailed anal- 
ysis of the various backgrounds shows that the most favourable 
energy range for detecting 7-rays from SNRs is probably from 
1 to 10 TeV. In addition, to characterise the spectra, detections 
in the GeV region and detections or upper bounds in the PeV 
region would be very useful and appear possible with current 
detectors. 

Galactic SNRs should be detectable 7-ray sources at dis- 
tances of up to a few kpc in regions of the ISM where the 

mean density is of order 1 cm~^ or more, if, as is generally 
supposed, they are the main sources of the cosmic rays below 
10" to lO^-"' eV. Even the detection of one or two SNR which 
are not dominated by a pulsar, at the level suggested in this 
paper would be a convincing demonstration of the SNR origin 
of cosmic rays. 
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Appendix A: The emissivity of very high energy neu- 
tral particles in p-p interactions 

On the assumption that the primary protons have a power- 
law differential energy spectrum, 

Np{E) = NoE-''+\ Al 

the emissivity for the secondary particles produced in p-p in- 
teractions may be written in the form 

£s{E) = {mx)%appcNp{E) A2 

where Opp is the inelastic p-p cross-section, c is the velocity of 
light and 

(ma;)g = / x"g{x) dx A3 
Jo 

is the so-called spectrum-weighted moment of the inclusive 
cross-section (see, e.g. Gaisser 1990; x is the ratio of the sec- 
ondary particle energy to that of the primary proton and ^(a;) 
is the dimensionless inclusive cross-section). 

A comprehensive discussion of the spectrum-weighted mo- 
ments for secondary hadrons based on the interaction of ac- 
celerator beams with fixed targets at beam energies < 1 TeV 
has been presented by Gaisser (1990). In table Al we present 
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the spectrum-weighted moments for TT^-mesons from Gaisser's 
book. 

As the production and decay of secondary tt" mesons is the 
main source of 7-rays (the contribution from the next most im- 
portant channel via rj mesons is only about {mx}^ / (mx)^o 
{m^o/rrir,)^ ^ 0.06) 



a + 1 



and thus 



A4 



gives for the first decay 



A9 



where 77 = 1 — m^/m^. An exact calculation of the i'^ spec- 
trum from decay (2) is more complicated, however it can be 

estimated as follows. In charged pion decays on average a frac- 
tion [1 — rj) of the pion energy is transferred to the muons, 
while approximately the same (~ 1/3) part of the muon en- 
ergy goes to the decay products, e*, V/j, and i/e- Thus in this 
simple approximation 



(ma;)" « {mx)"o 



Though the values of (mx)^o given in Table Al are cal- 
culated for fixed target accelerator data with E < 1 TeV it is 
expected that they correctly characterise also the energy region 

beyond 1 TeV (Gaisser, 1990). Indeed, the spectrum- weighted 
moments for gamma rays calculated on the basis of data from 
p-p colliders at = 630 GcV (the CERN UA7 experiment) 
can be approximated in simple form (A. Erlykin, personal com- 
munication) by 



log {rnx}" = 1.49 - 2.73a + 0.53a^ 



A6 



Comparison of Eq. A6 with (mx)'^ « (mx)JJo ^j^rj and the val- 
ues of (mx)^ in Table Al shows agreement within 20%, which 
seems very reasonable if account is taken of the uncertainty as- 
sociated with the contribution through »7-mesons, especially at 
high energies. 



A5 £^l>/£,^{a + l)e 



■■ 0.026 



AlO 



Equations AlO and A9 imply that both channels (1) and (2) 
make comparable contributions to the Vfj, production rate and 
that for hard proton spectra (ami) the ratio of gamma to 
neutrino production is close to unity. 

Over distances of order (_E/10^'^ eV) kpc neutrons can also 
be considered as quasistable neutral particles of possible astro- 
physical interest. The emissivity of neutrons can be calculated 
using the mean neutron multiplicity (m) « 1/4 and the nor- 



malised x-distribution of neutrons, g{x) = 3(1 - 
by Jones (1990) with the result 



, suggested 



(mx): 



3 r(3)r(Q-n) _ 3 1 

4 r(a-h3) ~ 2 (a-F3)(a-|-2)(a-|-l) 



All 



The various spectrum-weighted moments for stable neutral 
particles are presented in Table Al. 



Table Al. Spectrum weighted moments for the production of 
secondaries in p-p collisions. 



a 


1.0 


1.2 


1.4 


1.6 


1.8 


{mx}; Eq. (6) 


0.19 


0.094 


0.051 


0.030 


0.019 


{mx)'^o 

Gaisser (1990) 


0.17 


0.092 


0.066 


0.048 


0.036 


{mx): Eq. (11) 
Gaisser (1990) 


0.062 
0.05 


0.050 
0.05 


0.041 
0.04 


0.034 
0.03 


0.029 
0.03 


Eq. (9, 10) 
Gaisser (1990) 


0.95 
0.98 


0.80 
0.86 


0.67 
0.77 


0.56 
0.66 


0.46 
0.58 



At p-p interactions, as a result of the decay of charged 
pions, neutrinos will also be produced. From the point of view 
of detectability the most interesting are the muonic neutrinos, 
Vij, and Vp, produced in the decays (1) tt^ — > /i* -|- V/j, and (2) 

^ + fp + Ve- Taking into account that 

£^0 « ^£^± A7 
and the relation (eg Stecker, 1979) 

2^? JE/2r, E a + 1 



Appendix B: The diffuse Galactic 7-ray emission at 
very high energies 

The Galactic 7-ray background consists of two components, tt" 
and IC, resulting from interactions of cosmic ray nucleons and 
electrons with the ISM. The flux of the first component can be 
presented in the form (Aharonian, 1991) 

jf = ^£{> E)n (Nu) = (mx); c7ppr?JcR(> E) (uh) B1 

where app is the inelastic p-p cross-section, (mx); is the spec- 
trum weighted moment for 7-rays (see Appendix A), Jcr(> E) 
is the integral energy spectrum of the Galactic cosmic rays, r] is 
the contribution of nuclei, mostly *He, in the cosmic rays and 
the ISM to 77° production (for standard composition rj ^ 1.5, 
see eg Dermer, 1986) and (A^h) = (Nm) + (A^Hs) is the total 
hydrogen (atomic and molecular) column density in the given 
direction. 

Radio measurements in the 21 cm and 2.6 mm lines provide 
rather accurate estimates of (Nh), so the main uncertainty is 
associated with the cosmic ray flux above 1 TeV. Indeed, if the 
CR proton spectrum is taken to be 

Jcr(> £) « 7 X 10-"^ (p§v) ' cm-'s-^sr-i B2 
then at S > 100 MeV we have 

J^(> E) « 1.2 X 10-" (Nh)^^ cm-^s-^sr-i B3 
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whereas for the Hillas (1981) representation of the spectrum, 

E 



JcaO E) « 10- (^j 



-2 -1 -1 

cm s sr 



the corresponding diffuse 7-ray spectrum is 

E ^-1-6 



J^(> E) « 2.2 X 10" 



ITcV 



{Nh)22 cm ^s ^sr ^ B5 



where {Nh)22 = {^h) /lO^^ cm~^. Thus there is about a factor 
two uncertainty in the 7-ray flux above 1 TeV associated with 
uncertainty in the Galactic CR flux. 

The flux related to the electron component of the CR 
is even less certain. However the observations reported by 
Nishimura et al (1980, 1990) do extend to energies slightly 
above 1 TeV and show that the flux of high energy electrons is 
reasonably represented by a power-law. 



Je(S) = 7.6 X 10" 



E 



iTcv; 



cm ^s ^sr ^TeV 



B6 



in the region 10 GcV < E < 1 TeV. Here we will assume that 
this spectrum extends to well above 1 TeV. 

Although at 100 MeV the contributions from bremsstrahl- 
ung and IC processes on starlight and 2.7 K cosmic microwave 
background (CMB) are comparable, at VHE {E >! TeV) only 
IC on the CMB remains important because for 7-ray energies 
E < Arriic^ /ikT ^ 10'^ TeV it works in the Thomson limit 
and produces a hard spectrum of 7-rays. The emissivity of this 
process in the case of a narrow-spectrum background photon 
field with an average photon energy wo and number density no 
is (see, eg, Blumenthal and Gould, 1970) 



f-=.rgcnoiVo2^+3^!+M+^ 



(/3 + 3)(/3 + l)(/3 + 5) 



B7 



E 

meC^ 



where 



47r 



iVe(S) = -j-Je{E) = NoE- 



B8 



is the spectrum of the rclatlA'istic electrons. 

Thus, for IC of the high energy cosmic ray electrons on the 
CMB with Wo w 3kT and no w 400 cm~^ we have 



'^CMb(> E) — -— fcMB(> E)Lc 



4 X 10 



— V 

TeV; 



lOkpc 



B9 



-2 -1 -1 
cm s sr 



Equation B7 shows that, in the Thomson limit. 



B4 S-y oc noWo ^ WdWq ' 



BIO 



where Wo = nowo is the energy density of the photon field. 
Let us consider in this limit the IC emissivities for scattering 
the CMB, optical and far IR radiation fields. The optical ra- 
diation, with a mean photon energy of 1.5 eV is estimated to 
have an energy density of 0.5 eV cm~* and the far IR with a 
mean photon energy of 0.01 eV to have an energy density of 
0.2eVcm~^. Then in the Thomson limit 



tCMB 



0.15£:o^, 



c-lC 

tCMB 



Bll 



so that the 7-ray background due to the IC scattering of optical 
photons is roughly three times that due to the scattering of far 
IR photons which in turn is about twice that due to scattering 
of the CMB photons. However, the parameter S = E-^wo/rnic^ 
is « 6(£;T/lTeV) for optical photons and 0.04(£'t./1 TeV) 
for IR photons. The Thomson cross-section, ctt, is only ap- 
propriate for (5^1 and for 5 > O.Ol the Klein-Nishina cross- 
section, (tkn, must be used. Asymptotically 



0"KN 



1 - 2(5^/2 if 5 <C 1 

1^-1/2 ^1 + In (2J1/2)] if5>l 
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so that the cross-section drops by a factor of two at 5 ~ 0.05 
and by more than a factor 10 at 5 ~ 5. Thus at 1 TeV and 
above the dominant contribution to the IC background is in 
fact the component produced by the CMB (the spectra also 
steepen in the Klein-Nishina region) . 

It is also easy to show that the electron bremsstrahlung 
contribution is negligible when compared to the IC and 7r° 
backgrounds. 
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Finally, it is interesting to note that for reasonable confine- 
ment regions, Lc ~ lOkpc, the ratio of the fiux of IC 7-rays 
from the CMB to the primary electron flux is 



J^^(> E) 

Je(> E) 



0.012 
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— V 

TeVy 
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which is less than unity for E < 40 TeV. In fact this ratio 
is only an upper limit because Klein-Nishina effects become 
important for IC on the CMB at about 40 TeV (see above). 



where Le is the characteristic size of the confinement region of 
CR electrons in the Galaxy (in the given direction). 



This article was processed by the author using Springer- Verlag l^jX 
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